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OO ' ABSTRACT 

i We present the computation of effective refractive coefficients for in- 

homogeneous two-component grains with 3 kinds of inclusions with 
; minci = 3.0-h4.0i,2.0-h 1.0i,2.5-H0.0001i and a matrix with m^ = 1.33 O.Oli 

, for 11 volume fractions of inclusions from 0% to 50% and wavelengths A=0.5, 1.0, 2.0 

' and 5.0 /xm. The coefficients of extinction for these grains have been computed using 

I a discrete dipole approximation (DDA). Computation of the extinction by the same 

. method for grains composed of a matrix material with randomly embedded inclusions 

' has been carried out for different volume fractions of inclusions. A comparison of 

■<!;;j- \ extinction coefficients obtained for both models of grain materials allows to choose 

■ the best mixing rule for a mixture. In cases of inclusions with minci=2. 0-1-1. Oi and 

2.5-|-0.0001i the best fit for the whole wavelengths range and volume fractions of 
inclusions from to 50% has been obtained for Lichtenecker mixing rule. In case 
of minci = 3.0 -I- 4.0i the fit for the whole wavelength range and volume fractions of 
O \ inclusions from to 50% is not very significant but the best has been obtained for 

Hanai rule. For volume fractions of inclusion from to 15% a very good fit has been 
^ ■ obtained for the whole wavelength range for Rayleigh and Maxwell-Garnett mixing 



rules. 

Key words: ISM: general, dust, extinction, interstellar grains, mixing rules, discrete 
dipole approximation 



1 INTRODUCTION 

The primary goal of the research we present here has been 
to use computational electrodynamics methods to choose 
the best mixing rule for different materials (dielectric, semi- 
conductor and metal) of inclusions in dielectric matrix. 
Light scattering computation for composite particles can 
be made using the discrete dipole approximation (DDA). 
IVaidva et al.l ("2001) applied this method to calculate extinc- 
tion and scattering efficiencies for silicate sphere with em- 
bedd ed graphite inclusions a nd evaluated interstellar extinc- 
tionJ^idersei^^L| j|200^ tested the programs D DSCAT 
bv iDraine fc Flataul J2000ald) and M arCoDES bv iMarkell 
l|l998tl . The program by 'Markcl' (199?) is much faster than 
the one by iDraine fc Flatau (,2000a. .h) which is very time 
consuming, however, it is ment to be used mainly for cal- 
culating the extinction by sparse clusters, whereas for com- 
pact clusters the accuracy of this method is very low. For 
this reason the DDA methods are not used for mass study 
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of inter- or circumstellar extinction. This is why, many re- 
searchers still consider various models of grains such as host 
materi al in which other mater ials are embedde d (for ex- 
ample [Ma^i^^h^^ Jl989l) or iMaronl (119881) ') or core- 
mantle ("Jones? 198^) a nd mu ltilayer spherical particles 
( Voshchinnikov & Mathi^ il999ft ). Besides, the use of a par- 
ticular model compared to observatio ns can provide informa- 
tion about the inhomogeneous grains. IPerrin fc Lamvllll990l) 
computed cross section of extinction using: 1° DDA method 
assuming cubic lattice with elementary cell treated as inclu- 
sion of 10 A for spherical grain of radius 100 A and 2° Mie 
theory for grain with the same radius. They considered two 
mixtures: first a mat rix of si l icate with complex index of re- 
fraction proposed bv lDraiii3 il985ft and inclusions of thaolin 
(organic dielectric with imaginary part of refractive index 
k < 0.001) and the second case adesite containing inclusions 
of water ice. Both cases concern only a mixture of dielectrics 
for two values of volume fraction of inclusions in the infrared 
region. In their paper the effective optical constants applied 
to the Mie theory were obtained from the most popular 
rules in astrophysics: Maxwell-Garnett's and Bruggeman's. 
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IPerrin fc Lamvl lll99(J) compared cross sections of extinction 
obtained by both methods and stated that the apphcation 
of effective-medium theories (EMT) to the problem of hght 
scattering by inho mogeneous grains is not straightforward. 
Other authors fe.g. lChvlek et all teood) ) compared the DDA 
with EMT method but only for very limited range of volume 
fractions of inclusions, wavelengths and materials. 



2 MIXING RULES AND EFFECTIVE 
PERMITTIVITY 

We limited our study to the grains without electric charge, 
magnetic suspectibility and considered only two component 
mixtures where inclusions were monodisperse. 

In our investigation we tested six mixing rules, not only 
the most popular in astrophysics such as Maxwell-Garnett's 
and Bruggeman's. The considered mixing rules are presented 
below. 

The Maxwell-Garnett's l)Maxwell-Garnetlll9oll mixing 
rule, originally, has been applied to metal particles encap- 
suled in an insulating matrix and for small filling factors, 
now it is used without restriction. The f ollowing equation is 
known as the Maxwell-Garnett formula jBohren fc HuffmanI 
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For composite media consisting of spherical particles of both 
basic comp onents without the requirement of complete en- 
has established the following 
equation (Bohrcn & Ifuffmanl ll983 
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The Bruggeman's equation is symmetrical, i.e. Si and e„ 
can be exchanged. Simil ar equation, but not for spherical 



inclusions was derived bv lLandaueJ ^^5^) for conductivity 
of metallic mixtures. 

For spherical inclusions and random oriented ellipsoids 
Landau_fc Lifezi3 il960l) and later, along a different way, 
Loov^ngal ^1965^ have found the following relation: 
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The Hanai-Brugge man equation was originally derived by 
iBruggemanI (^33) for composite dielectric consisting of a 
host material with spherical inclusions. Hanai in 1961 mod- 
ified the Bruggeman's equation and generalised it for com- 
plex dielectric permittivities. The Hanai-Bruggeman equa- 
tion is following (eq.3.58 in Beek. I.1967.'l l: 
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Other mixing law for homogeneous and isotropic systems 
has been obtained bv iLichteneckeil (Il92(t l. The logarithmic 
mixing rule for a mixture of two components has the form 
of: 

loge = /logEi -f (1 - /)log£„. (5) 

iMeredith fc Tobiad (|l96ff l reconsidered Rayleigh's deriva- 
tion for a cubical array of spheres Ei enclosed in a medium 



£m and obtained an equation that appears to be more satis- 
factory at high values of /. The fo llowing equ ation describes 
improved Rayleigh's mixing rule (jBeeklll?- 
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In all formulas / is the volume fraction of inclusions and 
£m, Ei - the complex dielectric permittivity of a matrix and 
inclusion, respectively. The e (without a subscript) is the 
complex dielectric permittivity of the mixture. Generally the 
applicability of all mixing rules is restricted to low concen- 
tration of inclusions. In this paper we have extended the 
search for the best mixing rule up to 50% of concentration 
of inclusions. 



3 DETAILS OF CALCULATIONS 

We examined three different cases of the refractive index for 
inclusions rm^d = 3.0-|-4.0i, 2.0-l-l.Oj, 2.5+O.OOOli. The cho- 
sen refractive indices of inclusions refer t o those of materials 
examined bv lDraine fc GoodmanI il993l) . The matrix in all 
cases was dielectric with refractive index rrim ~ 1.33 + O.Oli. 
Because the refractive indices are not additive quantities 
they have been changed into complex permittivity relative 
to free space using the following formulae: 



(7) 



and 

e" = 2nk, (8) 

where: n and k are real and imaginary parts of complex re- 
fractive index, e' and e" are real and imaginary parts of com- 
plex permittivity. Next, we have calculated the complex per- 
mittivities for the mixture based on the above mixing rules. 
The computation has been carried out for volume fractions 
in the range from 0% to 50% with 5% step. The obtained 
complex permittivities have been changed into complex re- 
fractive indices using the following relations: 
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The extinction coefficients have been computed us- 
ing the Discrete Dipole Appr oximation, first proposed by 
iPurcell fc Pennvpackerl il973fl . DDA methods divide the 
particle into numerous polarizable elements of volume. The 
induced dipole polarizations in these cubes are determined 
self-consistently, then properties such as the extinction cross 
section are deter mined in terms o f the induced polariza- 
tion. In our studv. lPraine fc Flataul lll99'd) approach is used 
to compute the efficiency factor for extinction. The refrac- 
tion indices for each kind of mixture and all considered 
mixing rules were calculated using the FORTRAN routine 
DDSCAT.5alO given by Drainc fc Flatau (2000a). 

The efficiency factors for extinction Qf""™"^ for ho- 
mogeneous spherical grains (each element of DDA ar- 
ray has the same refractive index calculated for the mix- 
ture) have been computed for the following wavelengths: 
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A = 0.5, 1.0, 2.0, 5.0 nm and grains with radius r — 0.15 /im 
containing 1791 dipoles. The location of dipoles in the 
array for such pseudospherical particle has been ob- 
tai ned using the routine c alltarqet.f also made available 
bv iDraine fc Flataul (|2000a|). The subscript p in the sym- 
bol Q}^""^"^ corresponds to the appropriate mixing rule, j 
- the volume fraction of inclusion and the subscript 1 cor- 
responds to wavelength. The determined efficiency factors 
for extinction have been compared to those computed with 
the exact Mie theory. It has been found that the differ- 
ences betw een the t wo m ethods are very small similar to 
IPraine ~ GoodmanI (Il99,'^l results. 

We have applied a mixture of DDA elements with dif- 
ferent refractive indices while the ratio of number of DDA 
elements with refractive index of inclusion to all elements 
in the grain equals the volume fraction of inclusion. The 
assumed size of inclusion equals the volume of the DDA 
element (the size of inclusion cannot be smaller t han the 
elemen t volume considered in DDA approximation') . iDraind 
l|l985fl considered how surface granularity affects the accu- 
ra cy of calculatio ns. In our study validity criteria considered 
bv lDraiii3 il985l) are satisfied and furthermore the influence 
of surface granularity is reduced because it is almost the 
same in both cases: for elements with mixture of materi- 
als and for random mixture of elements with two kinds of 
materials. Similarly, the departure from sphericity for both 
pseudospheres is the same and does not influence the differ- 
ence of extinction in both cases. The arrangement of inclu- 
sions (DDA elements) is random. In order to obtain random 
number of DDA elements in the discrete dipole arrays we 
have used random number generator "Research Random- 
izer" available at http: //www. randomizer, org We have gen- 
erated 10 series of numbers corresponding to the given vol- 
ume fractions of inclusions out of all 1791 dipoles. The gen- 
erated numbers have been sorted in ascending order. The 
location of dipoles in the array for the spherical particle 
obtained from the routine calltarget.f was the same as for 
homogeneous grains with the only difference that for the 
randomly generated numbers the DDA elements had the 
refractive index of inclusions and for the others the refrac- 
tive index of a matrix. Figure Q presents the dependence of 
efficiency factors for extinction on the volume fraction for 
inclusions with regular and random locations. One can see 
from this figure that the geometrical structure of inhomo- 
geneous grain plays a very important role. Authors of the 
considered mixing rules have assumed a statistical distribu- 
tion of inclusions. Therefore, in our calculations the random 
distributions have been used. Because there is a scattering 
of results for different random distributions the efficiency 
factors for extinction have been calculated for 10 different 
distributions and then averaged. The dependence curves of 
extinction on volume fractions of inclusions have been in- 
terpolated by 5th degree polynomials. For grains with radii 
of r = O.lbfim 10 values of efficiency factors for extinction 
Qija'* depending on random location have been calculated 
for each of 4 wavelengths A=0.5, 1.0, 2.0 and 5.0 fim by us- 
ing the computer program DDSCAT.5alO. The subscript i 
in the symbol Qij""* corresponds to the number of random 
location of inclusions, j - the volume fraction of inclusion and 
the subscript 1 corresponds to wavelength. Next the mean 
extinction has been calculated as: 
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4 RESULTS AND DISCUSSION 

The figures I2I13I show the dependence of mean values of 
efficiency factors for extinction Qf^"'^ and Qf""™"^ on the 
volume fractions of inclusions for the given mixing rule and 
wavelength and the best fitted curves. 

The best fit of extinction for homogeneous grains con- 
sisting of material of averaged refractive index calculated 
from the mixing rule and extinction of grains consisting of 
DDA elements of two materials (matrix and inclusions) with 
different refractive indices (random) was obtained in the fol- 
lowing way: 

a) for the given mixing rule (subscript p) and the given 
wavelength (subscript 1) the Xp,; was calculated from 
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where summing is done for volume fractions of inclusions 
j, Q\J"^ is avera ged extinction coefficient for randomly lo- 
cated inclusions, is extinction coefficient for homo- 
geneous grains for a given mixing rule. Because the values 
of Q[""'' for different random locations of inclusions (i) was 
not treated as equivalent of measured values, therefore o"!j,p 
is not the standard deviation. (T!j,p is the fraction a of the 
averaged extinction coefficient Q\°'"'^. 
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The calculations of X;,p have been made for a =0.005, 0.01, 
0.025, 0.05 and 0.1. The distribution of Q\'^^f values has 
been assumed to be normal, 

b) The goodness of fit Aj.p has been calculated using the in- 
complete gamma fun ction Ai_p = gammq(0.5i', 0.5xf,p) de- 
scribed in jPress et al] (l993 ). where = N — M is the num- 
ber of degrees of freedom, N is the nuber of points in the 
curve (number of volume fractions of inclusions), M are ad- 
justable parameters (M=l number of mixing rules for calcu- 
lating the x?,p- The obtained values of fitting coefficients A;_p 
for each mixing rule and each wavelength for a = 0.025 are 
presented in Tables Q and |21 and for the metallic inclusions 
- a = 0.1 in Tables 13 and 0] because for smaller values of a 
the goodness-of-fit coefficients were very small. For the mix- 
ture of given materials determined by a mixing rule different 
fitting coefficients have been obtained for each wavelength. 
For different wavelength ranges different mixing rules may 
be used based on the given values of A;_p. 

c) In order to choose the best mixing rule in the whole con- 
sidered range of wavelengths the value 
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has been calculated and the procedure described in b) car- 
ried out leading to the obtained goodness-of-fit coefficient 
Ap. 

Summing up, we have examined 3 different materials as 
inclusions in a matrix of m^ ~ 1-33 + O.Oli: 
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Figure 1. Example of different values of extinction coefficients for different random locations of inclusions with respect to volume 
fractions of inclusions (dots and solid line best fit) and dependence of extinction on volume fraction for regular location of inclusions 
(stars and dashed line) both calculated using the DDA. Refractive index of inclusion is mine = 3.0 + 4.0i 



Table 1. Goodness of fit for each mixing rule and each wavelength A; p - mine = 3.0 + 4.0i, 
a = 10% of Q™d 



Mixing rule 


A;_p(0.5/im) 




A,^p(2.0Atm) 


A;_p(5.0^tm) 


Ap 


Bruggeman 


0.616E+00 


0.643E-40 


O.OOOE+00 


O.OOOE+00 


O.OOOE+00 


Hanai 


0.389E+00 


0.162E-06 


0.911E-02 


0.449E-01 


0.231E-06 


Lichtenecker 


0.996E+00 


0.750E-22 


O.OOOE+00 


O.OOOE+00 


O.OOOE+00 


Looyenga 


0.121E+00 


O.OOOE+00 


O.OOOE+00 


O.OOOE+00 


O.OOOE+00 


Maxwell— Garnctt 


0.240E-01 


0.595E-01 


0.535E-21 


0.616E-29 


O.OOOE+00 


Rayleigh 


0.132E+00 


0.433E-01 


0.120E-11 


0.123E-18 


0.129E-26 



(1) mine = 3.0 + 4. Oi (Figures l2l5ll 

(a) The values of goodness-of-fit coefficients for 
a =0.025 and 0.05 are too small to be considered for anal- 
ysis of the fit, therefore Tabled shows the values of Ai,p 
and Ap for a =0.1; 

(b) The goodness-of-fit of the extinction curve for each 
mixing rule strongly depends on the wavelength. For 
A = 0.5/im the best fit has been obtained for the Licht- 
enecker rule. However, for the whole range of wavelengths 
the Hanai rule gives the best fit for this material. It is nec- 
essary to point out that the value of the fitting coefficient 
Ap is very small. It may be necessary to investigate other 
mixing rules in order to obtain a better fit; 

(c) For a high volume factor of inclusions they may be 
placed so close that, on one hand, they may create inclu- 
sions of much bigger sizes and therefore more vulnerable 
to skin effect, on the other hand they may exceed the per- 
colation threshold. Therefore, we have considered the fit 
for the volume factor of inclusions in the range from to 
15 %. The results are presented in Table |5| from which we 



point out that for the whole wavelength range the best 
fit is obtained for Rayleigh and Maxwell-Garnett mixing 
rules. 



(2) mine = 2.0 + l.Oi (Figures I6l9t 

(a) The smallest value of a which gives the goodness- 
of-fit coefficients suitable for further analysis is 0.025 and 
for this value the coefficients Ai,p and Ap have been cal- 
culated; 

(b) For A = 0.5/im the best fit has been obtained for 
the Looyenga mixing rule (Fig|SJ. However, for the whole 
range of considered wavelengths the Lichtenecker rule 
gives the best fit (Table I^J. 

(3) mine = 2.5 + O.OOOli (Figures [THIT^ 

(a) The smallest value of a which gives the goodness- 
of-fit coefficients suitable for further analysis is 0.025 and 
for this value the coefficients Ai,p and Ap have been cal- 
culated; 

(b) For A = 0.5fim the best fit has been obtained for 
the Rayleigh and Maxwell-Garnett mixing rules. However, 
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Table 2. Goodness of fit for each mixing rule and each wavelength A; p - 
volume fraction of inclusions 0-15%, a = 10% of Q!!|?'* 



: 3.0 + 4.0i, 



Mixing rule 



A;_p(0.5/im) Ai.p(1.0/^m) Ai^p{2.0fj,m) A;_p(5.0/jm) 



Bruggeman 
Hanai 

Lichtenecker 
Looyenga 
Maxwell— Garnett 
Rayleigh 



0.826E+00 
0.592E+00 
0.999E+00 
0.157E-09 
O.lOOE+00 
O.lOOE+00 



0.383E-12 
0.182E-02 
0.158E-37 
O.OOOE+00 
0.997E+00 
0.997E+00 



0.518E-34 
0.184E-05 
O.OOOE+00 
O.OOOE+00 
0.599E+00 
0.624E+00 



O.OOOE+00 
0.597E-09 
O.OOOE+00 
O.OOOE+00 
0.399E+00 
0.423E+00 



O.OOOE+00 
0.950E-13 
O.OOOE+00 
O.OOOE+00 
0.994E+00 
0.996E+00 



Table 3. Goodness of fit for each mixing rule and each wavelength A; p - mine = 2.0 + l.Oi, 
cr = 2.5% of QS?'* 



Mixing rule 


Ai^piO.b/im) 


A,,p(1.0Atm) 


A,,p(2.0Atm) 


A;^p(5.0/im) 


Ap 




Bruggeman 


0.929E+00 


0.245E- 


-26 


O.OOOE+00 


O.OOOE+00 


O.OOOEH 


-00 


Hanai 


0.643E+00 


0.216E- 


-15 


0.420E-38 


O.OOOE+00 


O.OOOEH 


-00 


Lichtenecker 


0.256E+00 


0.237E4 


-00 


0.987E+00 


0.999E+00 


0.896EH 


-00 


Looyenga 


0.988E+00 


0.000E4 


-00 


O.OOOE+00 


O.OOOE+00 


O.OOOE+OO 


Maxwell— Garnett 


0.693E-01 


0.103E- 


-05 


0.946E-11 


0.450E-13 


0.317E- 


-26 


Rayleigh 


0.172E+00 


0.158E- 


-07 


0.131E-14 


0.327E-17 


0.736E- 


-35 



for the whole range of considered wavelengths the Lichte- 
necker rule gives the best fit (Table 3J. 

Considering the above results it is possible that differ- 
ent mixing rules should be applied for the same mixture of 
materials for different wavelengths, especially with higher 
volume fractions of inclusions. In spite of their deficiencies 
like simplifying idealisations while deriving them or basing 
only upon experiments, the mixing rules are able to provide 
important information about the inhomogeneous materials. 
In astrophysics generally the models by Maxwell-Garnett 
and Bruggeman are widely used because they are justified 
by theory. They are based on different topology of inclusions. 
Other mixing rules, less theoretically justified, are often ne- 
glected although they show better agr eement with experi- 
mental data. However, Zakr i et alj ()l998"l on the basis of ef- 
fective medium theory found physical grounds for the Licht- 
enecker rule which, as shown by the calculations, gives bet- 
ter fits for inclusions with refractive indices mine ~ 2.0 + l.Oi 
and mine = 2.5 + O.OOOli in the whole considered wavelength 
range and for mine = 3.0 + 4.0i in the short wavelength re- 
gion. 
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Table 4. Goodness of fit for each mixing rule and each wavelength A; p - 
mine = 2.5 + O.OOOli, a = 2.5% of Q™"* 



Mixing rule 


A;,p(0.5/im) 


Ai^p(1.0/^m) 




A;_p{5.0/jm) 


Ap 


Bruggeman 


0.259E-06 


0.371E-19 


0.483E-03 


0.342E-05 


0.143E-29 


Hanai 


0.177E+00 


0.645E-04 


0.277E-01 


0.999E+00 


0.138E-02 


Lichtenecker 


0.682E+00 


0.991E+00 


O.lOOE+01 


0.367E-01 


0.884E+00 


Looyenga 


0.217E-20 


0.455E-32 


O.llOE-03 


0.966E-12 


O.OOOE+00 


Maxwell— Garnett 


0.995E+00 


0.822E+00 


0.307E+00 


0.263E-10 


0.780E-05 


Rayleigh 


0.998E+00 


0.375E+00 


0.132E+00 


0.397E-03 


0.240E-01 
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Figure 2. The dependence of efficiency factors for extinction on 
the volume fraction of inclusions for different mixing rules and 
wavelengths. Refractive index of the matrix is 1.33 + O.Oli and 
for inclusion mine = 3.0 + 4.0i, the wavelength A = 0.5/im. The 
solid curve "random" shows the best fit for 10 extinction depen- 
dencies for randomly distributed DDA elements with refractive 
index corresponding to inclusion among elements with refractive 
index of a matrix. The grain radius in all cases is 0.15/.tm. 



Figure 3. Same as Fig.l^but for A = l.O^m. 
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Figure 4. Same as Fig. |2]but for A = 2.0/^m. 
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Figure 5. Same as Fig.|2]but for A = 5.0/^m. 
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Figure 6. Same as Fig. |2] but for mi, 
0.5/im. 



2.0 + l.Oi and A = Figure 7. Same as Fig. |2] but for mi, 

1.0/im. 
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Figure 8. Same as Fig. |2] but for mi, 
2.0^111. 



2.0 + l.Oi and A = Figure 9. Same as Fig. |2] but for mi, 

5.0/im. 



2.0 + l.Oi and A 
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Figure 10. Same as Fig.|2lbut for mine = 2.5 + O.OOOli and A = Figure 11. Same as Fig.|2lbut for mine = 2.5 + O.OOOli and A = 

0.5/im. l.Ofim. 
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Figure 12. Same as Fig. l2lbut for mi, 
2.0Mm. 



: 2.5 + O.OOOli and A : 



Figure 13. Same as Fig. l2lbut for mi, 
5.0/im. 
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